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An  unconventional  nonlinear  elastic  theory  is  advocated  for  solids  undergoing  large  com¬ 
pression  as  may  occur  in  shock  loading.  This  theory  incorporates  an  Eulerian  strain 
measure,  in  locally  unstressed  material  coordinates.  Analytical  predictions  of  this  theory 
and  conventional  Lagrangian  theory  for  elastic  shock  stress  in  anisotropic  single  crys¬ 
tals  of  aluminum,  copper  and  magnesium  are  compared.  Eulerian  solutions  demonstrate 
greater  accuracy  compared  to  atomic  simulation  (aluminum)  and  faster  convergence 
with  increasing  order  of  elastic  constants  entering  the  internal  energy.  A  thermomechan¬ 
ical  framework  incorporating  this  Eulerian  strain  and  accounting  for  elastic  and  plastic 
deformations  is  outlined  in  parallel  with  equations  for  Lagrangian  finite  strain  crys¬ 
tal  plasticity.  For  several  symmetric  crystal  orientations,  predicted  values  of  volumetric 
compression  at  the  Hugoniot  elastic  limit  of  the  two  theories  begin  to  differ  substan¬ 
tially  when  octahedral  or  prismatic  slip  system  strengths  exceed  about  1%  of  the  shear 
modulus.  Predicted  pressures  differ  substantially  for  volumetric  compression  in  excess  of 
5%.  Predictions  of  Eulerian  theory  are  closer  to  experimental  shock  data  for  aluminum, 
copper,  and  magnesium  polycrystals. 

Keywords:  Elasticity;  plasticity;  shock  physics;  crystals;  metals. 


1.  Introduction 

Nonlinear  elastic  constitutive  models  are  needed  if  strains  or  rotations  are  large 
in  anisotropic  materials,  or  if  the  stress-strain  response  of  the  material  is  suffi¬ 
ciently  nonlinear.  Thermodynamic  potentials  of  usual  Lagrangian  theories  of  non¬ 
linear  elasticity  of  crystals  [Wallace,  1972;  Thurston,  1974;  Teodosiu,  1982;  Clayton, 
2011]  incorporate  the  right  Cauchy-Green  strain  tensor  E,  a  function  of  the  (elas¬ 
tic)  deformation  gradient  and  its  transpose.  Such  theories  have  been  used  with 
some  success  for  describing  high  rate  behavior  of  ceramics  and  minerals  [Winey  and 
Gupta,  2004;  Clayton,  2009,  2010;  Foulk  and  Vogler,  2010],  metals  [Clayton,  2005; 
Winey  and  Gupta,  2006;  Vogler  and  Clayton,  2008],  and  concrete  [Clayton,  2008], 
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to  compressive  strains  of  up  to  several  percent  or  more.  Another  recent  model  for 
dynamic  deformation  of  cubic  single  crystals  incorporates  the  deviatoric  part  of  the 
elastic  Green  strain  for  the  elastic  shear  response  coupled  with  a  pressure-volume- 
temperature  equation-of-state  (EOS)  [Luscher  et  al. ,  2013]  for  the  pressure.  For 
isotropic  (e.g.,  untextured  polycrystalline)  solids,  nonlinear  elasticity  incorporating 
the  logarithm  of  the  left  stretch  (Hencky’s  tensor)  has  been  found  accurate  for  a 
number  of  materials  subjected  to  moderately  large  strains  [Anand,  1979]. 

In  the  present  work,  alternative  theory  is  advanced  incorporating  an  uncon¬ 
ventional  Eulerian  strain  measure  D,  with  components  referred  to  locally  relaxed 
material  coordinates,  that  is  a  function  of  the  inverse  elastic  deformation  gradient 
and  its  transpose.  This  strain,  which  has  historically  shown  promise  for  hydrostatic 
compression  [Davies,  1974;  Weaver,  1976;  Perrin  and  Delannoy,  1978],  is  applied 
in  this  work  to  problems  involving  simultaneous  volume  change  and  shear,  i.e. , 
shock  compression.  Herein,  the  designation  “Eulerian”  is  used  to  refer  to  a  strain 
that  is  a  function  of  the  inverse  deformation  gradient  [Davies,  1974;  Weaver,  1976; 
Nielsen,  1986],  and  not  necessarily  one  referred  to  spatial  coordinates.  Strain  tensor 
D  [defined  formally  in  Eq.  (2.3)]  is  thus  Eulerian  but  referred  to  material  coordi¬ 
nates,  the  latter  feature  enabling  its  use  in  constitutive  models  of  anisotropic  solids. 
In  contrast,  the  Almansi  strain  tensor  e  entering  Murnaghan’s  theory  for  isotropic 
solids  [Murnaghan,  1937]  is  simultaneously  Eulerian  and  referred  to  spatial  coordi¬ 
nates,  as  is  Hencky’s  strain  measure  [Anand,  1979]. 

In  recent  work  [Clayton,  2013],  a  complete  thermoelastic  theory  incorporating 
D  was  developed  in  parallel  with  Lagrangian  E-based  theory,  and  analytical  solu¬ 
tions  to  the  one-dimensional  (planar)  shock  problem  were  derived.  Predictions  of 
these  solutions  were  compared  for  single  crystals  of  quartz,  sapphire,  and  diamond 
shocked  in  pure  mode  directions. a  These  specific  results  did  not  demonstrate  a  clear 
advantage  of  one  formulation  over  the  other,  and  elastic  constants  of  order  three 
(sapphire)  and  order  four  (quartz  and  diamond)  were  needed  to  match  Hugoniot 
data.  On  the  other  hand,  analytical  solutions  for  hydrostatic  compression,  uniaxial 
compression,  and  simple  shear  of  generic  cubic  crystals  with  Cauchy  symmetry  and 
a  characteristic  pressure  derivative  of  the  ambient  bulk  modulus  of  four  demon¬ 
strated  certain  advantages  of  Eulerian  theory  over  Lagrangian  theory  in  terms  of 
general  accuracy  and  intrinsic  stability  [Clayton,  2013]. 

In  Sec.  2.3,  this  thermoelastic  theory  is  newly  applied  to  analyze  shock  compres¬ 
sion  of  metals.  Since  the  experimental  Hugoniot  elastic  limit  (HEL)  usually  occurs 
at  small  uniaxial  compressive  strain  in  ductile  metals,  purely  elastic  predictions  of 
material  strength  that  omit  relaxation  due  to  slip  and/or  twinning  are  idealizations 
when  deformations  are  finite,  though  large  elastic  deformation  in  the  absence  of  slip 


aDirections  for  purely  longitudinal  shock  propagation  were  normal  to  X-  and  Z-cut  faces  in  quartz 
and  sapphire  and  normal  to  an  X-cut  face  diamond;  in  that  analysis  of  K-cut  quartz,  transverse 
particle  motion  was  neglected. 
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may  be  possible  for  small  volumes  of  material  or  near  defect  cores  [Teodosiu,  1982; 
Clayton  et  al. ,  2014].  Apart  from  purely  theoretical  interest,  application  of  analyti¬ 
cal  thermoelastic  solutions  to  metals  is  useful  for  comparison  with  or  validation  of 
molecular  dynamics  simulations  or  their  interatomic  potentials  [Zimmerman  et  al. , 
2011],  since  defect-free  configurations  can  be  shocked  to  finite  strains  in  simulations 
of  small  volumes  over  short  time  scales  in  such  controlled  simulations.  Results  are 
reported  here  for  aluminum,  copper  and  magnesium. 

Crystal  plasticity  theory  is  generally  needed  for  accurate  continuum  modeling  of 
strength  of  single  crystals  shocked  above  their  HEL.  In  Sec.  3.1,  a  thermomechan¬ 
ical  framework  accounting  for  finite  elastic  and  plastic  deformations  is  formulated, 
incorporating  Eulerian  elastic  strain  D.  For  purposes  of  comparison,  corresponding 
equations  for  Lagrangian  theory  incorporating  strain  tensor  E  are  also  reported. 
Analytical  results  are  available  for  isentropic  compression  of  symmetric  orientations 
of  single  crystals  described  by  linear,  rate-independent  crystal  plasticity  [Johnson, 
1972;  Johnson,  1974],  but  these  exclude  geometric  and  material  nonlinearities  and 
thermal  effects.  Finite  strain  analytical  solutions  are  not  available  for  crystals  with 
anisotropic  inelasticity,  and  analysis  of  nonlinear  elastic-plastic  shock  behavior  can 
be  complex  even  for  isotropic  solids  [Germain  and  Lee,  1973;  Wallace,  1980].  In 
Sec.  3.2,  predictions  of  yield  behavior  (resolved  shear  stresses  on  common  slip  sys¬ 
tems)  and  temperature  rise  under  shock  compression  are  compared  for  generic  solids 
with  properties  representative  of  metallic  crystals,  in  order  to  quantify  strain  regimes 
wherein  the  choice  of  nonlinear  elastic  model  becomes  important.  Predictions  of  the 
pressure-volume  EOS  obtained  from  Eulerian  and  Lagrangian  theories  are  com¬ 
pared  with  shock  data  [Marsh,  1980],  demonstrating  greater  accuracy  of  the  former 
at  high  pressures  if  material  strength  can  be  neglected. 


2.  Nonlinear  Elasticity 

Section  2.1  reviews  continuum  elasticity  theory  incorporating  Lagrangian  and  Eule¬ 
rian  strain  measures.  Section  2.2  presents  analytical  solutions  for  steady,  planar 
elastic  shock  propagation  obtained  from  each  theory.  Only  essential  aspects  are 
reported  here  since  complete  derivations  have  been  recently  presented  elsewhere 
[Clayton,  2013].  Section  2.3  newly  applies  these  solutions  to  single  crystalline  alu¬ 
minum,  copper,  and  magnesium. 


2.1.  Continuum  thermoelasticity 

Letting  t  denote  time,  spatial  and  initial  material  coordinates  are  related  by  the 
motion 


x  =  x(X,  t). 


(2.1) 
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The  deformation  gradient  F  and  volume  ratio  J  are 

F  =  <9x/<9X  =  V0x,  J  =  det  F  >  0.  (2.2) 

Lagrangian  (E)  and  Eulerian  (D)  strain  tensors,  each  with  components  referred  to 
a  Cartesian  reference  coordinate  system,  are  defined  by 

E=i(FTF-l),  D  =  ^(1  -  F-1F-r).  (2.3) 

Like  E,  tensor  D  is  invariant  under  spatial  rotation,  and  it  differs  from  Almansi 
strain  e  =  i(l  —  F-TF-1)  of  Murnaghan’s  theory  [Murnaghan,  1937]  that  is  not 
rotationally  invariant  and  is  thus  restricted  (with  regards  to  use  as  a  thermodynamic 
state  variable)  to  isotropic  elasticity.  In  the  smooth  case,  local  momentum  and 
energy  balances  are,  for  no  body  forces  or  heat  transfer 

V  ■  er  =  pi),  er  =  crT ,  U  =  Jer  :  Vv.  (2-4) 

Cauchy  stress  is  a,  internal  energy  per  unit  reference  volume  is  U,  particle  velocity 
is  U,  and  spatial  mass  density  is  p  =  po/J.  Denote  by  4'  =  U  —  dp  the  Helmholtz  free 
energy  density,  with  9  and  p  denoting  temperature  and  entropy.  Thermodynamic 
potentials  for  Lagrangian  and  Eulerian  theories  are,  respectively 

U  =  U{  E,  17),  U=U(  D,  77);  4>  =  4>(E,  9),  T  =  <L(D,  (9).  (2.5) 

For  description  of  anisotropic  elasticity  in  crystals,  strain  energy  potentials  associ¬ 
ated  with  (2.5)  are  typically  expanded  as  Taylor  polynomials  in  strain  tensors  E 
(Lagrangian)  or  D  (Eulerian),  with  the  order  of  the  polynomial  corresponding  to 
the  maximum  order  of  elastic  constant  present  in  the  representation  (e.g.,  a  third- 
order  polynomial  incorporates  elastic  constants  of  up  to  order  three) .  As  explained 
in  the  Appendix,  with  increasing  order  such  polynomials  should  ultimately  converge 
to  the  same  value  of  energy  at  a  given  strain  of  effective  magnitude  less  than  unity. 
Differences  between  predictions  of  Lagrangian  and  Eulerian  representations  in  (2.5) 
result  from  truncation  of  their  Taylor  polynomials  at  a  certain  order. 

Thermoelastic  constitutive  equalities  and  temperature  change  given  by  the  local 
energy  balance  are  derived  using  standard  principles  [Clayton,  2013].  For  Lagrangian 
and  Eulerian  theory,  these  are,  respectively 

S  =  dU/dF  =  d^/dF  =  JF^ctF-7’,  9  =  dU/dp, 

8S 

p  =  —chi'/ 89:  c9  =  9—:  E,  (2.6) 

S  =  dU/dV  =  di>/d  D  =  JFTaF ,  9  =  dU/dp , 

<9S 

p  =  -d^/d9-  c9  =  9—  :  D.  (2.7) 
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Specific  heats  per  unit  reference  volume  at  fixed  strain  are  c  =  dU/dd  and  c  = 
dU/dO. 

As  noted  by  Davies  [1973],  Eulerian  strain  tensor  D  was  perhaps  first  intro¬ 
duced  by  Thomsen  [1972]  for  studying  compression  of  cubic  crystals.  For  cubic  or 
isotropic  solids  subjected  to  spherical  deformation  of  the  form  F  =  J1/,31,  stress 
er  =  —pi  is  hydrostatic  with  p  the  Cauchy  pressure.  Under  these  (isentropic)  condi¬ 
tions,  as  outlined  in  the  Appendix,  Lagrangian  and  Eulerian  theories  with  internal 
energy  functions  (2.5)i  expanded  as  third-order  Taylor  polynomials  in  E  and  D, 
respectively,  yield  the  following  respective  pressure- volume  equations-of-state: 


p  =  - dU/dJ  =  -B0  (J"1/3  -  J1/3) 


1-^B'(J2/3-!) 


p  =  -dU/dJ  =  -B0  (J_7/3  -  J_5/3) 


1+^(B'  -4)(J-2/3-!) 


(2.8) 

(2.9) 


where  Bo  and  Bq  are  the  isentropic  bulk  modulus  and  its  pressure  derivative  in  the 
reference  state.  Equation  (2.9)  is  identical  to  the  two-parameter  Birch-Murnaghan 
EOS  [Birch,  1947;  Thomsen,  1970].  As  noted  in  the  Appendix,  differences  in  pres¬ 
sures  predicted  by  (2.8)  and  (2.9)  at  the  same  volume  ratio  J  arise  from  truncation 
of  the  associated  internal  energy  functions  at  order  three  in  either  finite  strain  mea¬ 
sure.  Such  differences  would  presumably  eventually  disappear  with  the  addition  of 
a  sufficient  number  of  higher-order  terms  involving  higher-order  derivatives  Bg  etc. 
[Birch,  1978], 


2.2.  Analysis  of  planar  shock  compression 

A  steady  planar  shock  with  natural  velocity  T>  moving  through  an  unstressed  solid 
can  be  described  by  the  Rankine-Hugoniot  equations  [Germain  and  Lee,  1973; 
Thurston,  1974] 

V  =  v/(1-J),  P  =  PoVv,  IU}  =  ^pQv2.  (2.10) 

Particle  velocity  is  v  and  axial  shock  stress  is  P,  positive  in  compression.  Material 
ahead  of  the  shock  is  assumed  at  rest,  and  the  jump  in  a  quantity  is  its  value  in  the 
shocked  state  minus  its  value  in  the  initial  state.  For  compression  (uniaxial  strain) 
along  X± 


P=-on,  J  =  Fn  =  (1  +  2P)1/2  =  (1  —  2D)-1/2.  (2.11) 

Consider  the  response  of  an  elastic  crystal,  without  defects.  When  the  internal 
energy  is  a  linear  function  of  entropy,  analytical  solutions  to  the  planar  shock  prob¬ 
lem  can  be  derived  [Clayton,  2013].  Internal  energy  functions  (to  order  four  in  strain) 
and  conjugate  thermodynamic  stresses  for  Lagrangian  and  Eulerian  theories,  for 
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uniaxial  strain,  are 

U  =  iCuE2  +  ^CniE3 

+  iCmi EA  -  90  (r^  +  if n E2  -ljv,  S=  -J^P  =  dU/8E, 

(2.12) 

V  =  icn£>2  +  icmD3 

+  ^-jCim-D4  —  0o  (tiD  +  —TiiD2 

(2.13) 


-  1 )  n,  S=  ~J3P  =  dU/dD , 


where  0  =  0o  >  0  and  U  =  rj  =  0  are  state  data  for  material  ahead  of  the  shock. 
Second-order  elastic  constant  Cn  and  first-order  Griineisen  parameter  Id  are  equal 
for  Eulerian  and  Lagrangian  theories.  Third-order  elastic  constants  and  second- 
order  Griineisen  parameters  are  related  by  [Weaver,  1976;  Perrin  and  Delannoy, 
1978;  Clayton,  2013] 

Cm  =  Cm  +  12Cn,  fii  =  fn  +  TTi.  (2-14) 


Such  relationships  among  elastic  constants  are  addressed  in  more  detail  in  the 
Appendix.  Simultaneous  solution  of  (2.10)-(2.13)  with  use  of  binomial  series  expan¬ 
sions  for  Fn(E)  and  Fn(D)  and  leads  to  the  following  fifth-order  polynomials  for 
entropy  r/  generated  across  the  shock  [Clayton,  2013] 

5  5 

0 (E)  =  J2  a^Ek,  ti(D)  =  J2  bkDk,  (2.15) 

fc= 0  fc= o 

ao  =  b0  =  ai  =  &i  =  a2  =  b2  =  0,  (2-16) 

°3  =  i20o  ^ 111  +  3("n^  ~  (2-17) 


a4  —  —  0O  1  [Cim  +  3Cm  —  6Cn  +  ri(Cm  +  3Cn)], 

1  _  A  „  (2-18) 
64  =  — 0O  [C1111  —  9Cm  —  6Cn  +  ri(Cm  —  9Cn)], 

<25  =  — 0O  1[2Cmi  —  6Cm  +  15Cn 

4o 

+  ri(Cmi  +  2Cm  —  9Cn)  +  r((Cm  +  3Cn)], 

(2.19) 

05  =  ^^0  1[—  6C1111  —  6Cm  —  9Cn 

+  ri(Cmi  —  6Cm  —  33Cn)  +  rf  (Cm  —  9Cn)]. 
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Lagrangian  solution  77(E)  was  derived  by  Thurston  [1974].  Eulerian  solution  77(E) 
has  apparently  been  derived  first  by  the  present  author  [Clayton,  2013].  Conjugate 
stresses  are,  to  order  five  in  strain  measure  E  or  D, 

S  =  Ch-E  +  -Cm-E2  +  ^-Cnn  —  0oria,3^  E3 

—  e0E4[(Tia4  +  Fna3)  +  flofTids  +  rna4)£],  (2.20) 

S  =  Cn D  +  —Cm D2  +  ^-Cini  —  6»0ri63^  D3 

-  0oC4[(r164  +  r  1163)  +  (ri&s  +  f  11  bi)D).  (2.21) 

Shock  stress  P  is  then  obtained  from  these  stresses  using  (2.11),  (2.12)2,  and  (2.13)2 
and  is  fully  known  as  a  function  of  compression  ratio  J  =  V/Vq.  For  Lagrangian 
theory 

P  =  —(1  +  2E)1^2  jcii£  +  -Cm-E2  +  ^-Cnii  —  0oFid3^  E3 

-  0o£4[(ria4  +  f na3)  +  0o(ria5  +  rua4)£]},  (2.22) 

where  E,  D  and  J  are  related  through  (2.11)2-  For  Eulerian  theory 

P  =  -(1  -  2C)3/2  |Cii£>  +  icm D2  +  Qcmi  -  D3 

-  0o^4[(ri64  +  f  1163)  +  (rifts  +  f  n64)£>]|.  (2.23) 

Shock  velocity  T>  and  particle  velocity  v  can  then  be  found  using  (2.10)i  and  (2.10)2- 


2.3.  Application  to  metallic  single  crystals 

Properties  for  aluminum  (Al),  copper  (Cu)  (both  with  cubic  crystal  structure, 
shocked  along  cube  axis  [100])  and  magnesium  (Mg)  (hexagonal  crystal  structure, 
shocked  normal  [a]  and  parallel  [c]  to  the  c-axis  [0001])  are  listed  in  Table  1.  These 
cases  all  represent  so-called  pure  mode  directions  wherein  planar  impact  with  no 
transmitted  shear  stress  results  in  a  plane  wave  with  a  purely  longitudinal  com¬ 
ponent,  so  the  one-dimensional  analysis  of  Sec.  2.2  applies  exactly.  Second-  and 
third-order  elastic  constants  for  Al  [Thomas,  1968],  Cu  [Hiki  and  Granato,  1966], 
and  Mg  [Slutsky  and  Garland,  1957;  Naimon,  1971]  are  isentropic.  Third-order  con¬ 
stants  are  converted  from  mixed  coefficients  measured  ultrasonically.  Lagrangian 
fourth-order  constants  are  obtained  from  atomic  theory  for  Cu  [Wang  and  Li,  2009] 
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Table  1.  Single  crystal  properties  (6q  =  295  K,  Cap  in  GPa 
and  po  in  g/cm3). 


Property 

Al  [100] 

Cu  [100] 

Mg  [a-axis] 

Mg  [c-axis] 

Cn 

107 

166 

59.4 

61.6 

Cm 

-1080 

-1279 

-664 

-728 

Cm 

203 

715 

49 

12 

Cnn 

25000 

11900 

8170 

7380 

Cnn 

10500 

2000 

1220 

893 

Id 

2.17 

1.97 

1.52 

1.52 

PO 

2.70 

8.96 

1.74 

1.74 

and  Mg  [Rao  and  Padmaja,  1990],  and  from  experimental  shock  data  for  A1  [Zim¬ 
merman  et  al.,  2011].  Griineisen  parameters  are  obtained  from  the  identity 

6 

ri  =  £  Cr/ja^/cp,  (2.24) 

3-1 

with  a/3  (Voigt  notation)  and  cp  thermal  expansion  coefficients  and  specific  heat  at 
constant  pressure  [Wallace,  1972].  Higher-order  Griineisen  parameter  fn  is  found 
assuming  pT p  =  po?o/3  =  constant  [Wallace,  1980],  leading  to  fn  =  Id  [Clayton, 
2013] .  Note  that  for  c-axis  compression  of  Mg,  notation  of  Table  1  implies  C33  — >  Cn 
is  the  longitudinal  stiffness  normal  to  the  basal  plane,  and  similarly  for  higher-order 
constants. 

For  each  crystal,  predictions  for  shock  stress  versus  volume  ratio  are  made  using 
the  analytical  solutions  of  Sec.  2.2  for  Eulerian  and  Lagrangian  theories,  (2.22)  and 
(2.23),  respectively.  Model  predictions  of  shock  stress  P  normalized  by  second-order 
isentropic  elastic  constant  Cn  are  shown  in  Fig.  1  for  Al,  Fig.  2  for  Cu,  Fig.  3  for 
Mg  [a-axis]  and  Fig.  4  for  Mg  [c-axis].  Elastic  constants  of  up  to  order  four  are 
considered  in  results  labeled  “4th  order” .  For  Al  and  Mg,  the  relationship  [Clayton, 
2013] 

C1111  =  Cun  —  18Cin  —  318Cn,  (2.25) 

obtained  (see  Appendix)  by  assuming  diU/dFf1  =  (FU/dF^  at  the  initial  reference 
state,  was  used  to  determine  the  Eulerian  fourth-order  constant.  This  assumption 
was  unsatisfactory  for  Cu,  so  for  that  material,  Cnn  was  matched  to  the  Lagrangian 
4th  order  solution.  Results  labeled  “3rd  order”  and  “2nd  order”  are  obtained,  respec¬ 
tively,  by  setting  fourth-order  and  both  third-  and  fourth-order  elastic  constants  to 
zero. 

All  longitudinal  higher-order  elastic  constants  (i.e. ,  all  third-  and  fourth-order 
constants)  are  smaller  in  magnitude  for  Eulerian  than  Lagrangian  theory  for  these 
metals,  as  is  evident  from  Table  1.  Furthermore,  shock  stress  predictions  of  2nd 
and  3rd  order  models  are  generally  closer  to  those  of  4th  order  theory  for  Eule¬ 
rian  anisotropic  elasticity  than  Lagrangian  anisotropic  elasticity,  as  is  evident  from 
Figs.  1-4.  Table  2  quantifies  error  (%)  of  2nd  and  3rd  order  predictions  rela- 
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1.00  0.98  0.96  0.94  0.92  0.90 


vi  vn 


Fig.  1.  Axial  shock  stress  in  aluminum  single  crystal  shocked  elastically  along  [100]  cube  axis; 
atomic  simulation  data  from  [Zimmerman  et  al.,  2011]. 


Fig.  2.  Axial  shock  stress  in  copper  single  crystal  shocked  elastically  along  [100]  cube  axis. 


tive  to  4th  order  predictions,  computed  as  2-(2nd  or  3rd  order  result  —4th  order 
result)/(2nd  or  3rd  order  result  +  4th  order  result).  For  each  crystal  type,  such 
errors  are  almost  always  smaller  in  magnitude  for  Eulerian  theory  (E2,  E3)  than 
for  Lagrangian  theory  (L2,  L3)  at  a  given  volume  ratio  and  order  of  approximation. 
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Fig.  3.  Axial  shock  stress  in  magnesium  single  crystal  shocked  elastically  along  [1210]  a-axis. 


Fig.  4.  Axial  shock  stress  in  magnesium  single  crystal  shocked  elastically  along  [0001]  c-axis. 


Together,  these  observations  imply  a  faster  converging  series  in  (2.23)  than  in  (2.22) 
as  the  number  of  higher-order  constants  is  increased,  and  greater  accuracy  of  Eule- 
rian  theory  than  Lagrangian  theory  when  the  same  number  of  elastic  constants  (i.e., 
the  same  order  of  Taylor  polynomial)  is  used  for  each  internal  energy  function  in 
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Table  2.  Relative  error  (%)  in  shock  stress  P  predicted  by  2nd  (L2,  E2)  and  3rd  (L3,  E3)  order 
theories. 


V/V0 

Aluminum  [100] 

Copper  [100] 

Magnesium  [0001]* 

L3 

E3 

L2 

E2 

L3 

E3 

L2 

E2 

L3 

E3 

L2  E2 

0.96 

-5.0 

-3.1 

-23.1 

+1.1 

-1.6 

-0.4 

-15.8 

+9.3 

-2.5 

-0.4 

-23.4  -0.0 

0.92 

-15.7 

-14.1 

-48.2 

-5.0 

-5.4 

-2.1 

-31.6 

+  19.9 

-7.9 

-2.0 

-45.1  -1.2 

0.88 

-28.1 

-34.5 

-71.2 

-19.8 

-10.4 

-6.2 

-46.5 

+31.5 

-14.5 

-5.2 

-64.2  -3.8 

*Similar  trends  (i.e.,  smaller  error  for  Eulerian  theory)  are  obtained  for  compression  along  a-axis 
(not  shown). 


(2.5)i.  For  Al,  as  shown  in  Fig.  1,  the  analytical  solution  incorporating  4th  order 
Eulerian  theory  also  more  closely  matches  atomic  predictions  [Zimmerman  et  al., 
2011]  of  elastic  shock  stress,  with  analogous  4th  order  Lagrangian  theory,  and  all 
lower  order  theories,  apparently  too  compliant.  At  J  =  0.923,  the  error  in  4th  order 
Lagrangian  theory  relative  to  atomic  simulation  is  —7.4%,  compared  to  only  —1.4% 
for  4th  order  Eulerian  theory. 

While  these  results  suggest  apparent  advantages  of  Eulerian  over  Lagrangian 
theory  for  shock  compression  of  metallic  single  crystals,  it  is  cautioned  that  the 
predictions  are  purely  (thermo)elastic  and  therefore  are  strictly  applicable  only 
for  very  small  volumes  —  such  as  in  the  aforementioned  atomic  simulations  or 
in  the  immediate  vicinity  of  pinned  defect  cores  —  wherein  plastic  deformation 
does  not  occur.  In  laboratory  scale  specimens,  yielding  would  commence  in  each  of 
these  pure  metals  at  small  compressions  at  which  effects  of  higher-order  constants 
and  differences  in  shock  stresses  predicted  by  the  two  theories  would  be  nearly 
negligible.  However,  nonlinear  elastic  effects  on  deviatoric  stresses  could  still  be 
important  at  larger  compressions  after  yielding,  particularly  for  lower  symmetry 
materials  with  restricted  numbers  of  slip  planes/directions  [Johnson,  1972],  and 
the  nonlinearity  in  pressure-volume  response  always  becomes  important  at  larger 
compressions  regardless  of  shear  strength,  as  will  be  demonstrated  explicitly  later. 


3.  Nonlinear  Elastic— Plastic  Theory 

Nonlinear  crystal  elastoplasticity  theory  is  needed  to  address  the  deviatoric  response 
of  anisotropic  crystals  when  shear  stresses  become  large  enough  for  defect  motion 
to  occur.  General  continuum  frameworks  are  outlined  in  Sec.  3.1  wherein  elastic¬ 
ity  is  described  by  Lagrangian  or  Eulerian  strain  measures,  and  where  plasticity  is 
described  by  shearing  on  preferred  planes  in  preferred  directions.  Thermodynamic 
results  are  somewhat  standard  for  crystal  plasticity  with  Lagrangian  elasticity  [Clay¬ 
ton,  2011],  but  are  new  for  the  present  kind  of  Eulerian  elasticity.  In  the  context  of 
shock  loading,  a  crystal  plasticity  framework  is  also  required  for  proper  derivation  of 
driving  forces  (i.e.,  conjugate  shear  stresses)  required  for  slip,  twinning,  or  possible 
shear  fracture  in  crystals  at  the  onset  of  yield  (i.e.,  at  the  HEL)  and  at  higher  shock 
stresses.  Analysis  of  such  phenomena  is  undertaken  in  Sec.  3.2  via  representative 
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calculations,  followed  by  further  consideration  of  the  volumetric  response  predicted 
by  Lagrangian  and  Eulerian  models. 

3.1.  Continuum  crystal  plasticity 

Deformation  gradient  F  of  (2.2)  is  decomposed  multiplicatively  into  elastic  and 
inelastic  (i.e.,  plastic)  parts 

F  =  FPFP,  Je  =  detFE  >  0,  Jp=detFp>0.  (3.1) 

In  general,  neither  FE  nor  Fp  is  integrable  to  a  continuous  motion  field,  i.e.,  these 
are  anholonomic  deformations  [Clayton,  2012].  The  time  rate  of  plastic  deformation 
is  dictated  by  the  usual  slip  kinematics  [Asaro,  1983] 

FPFP_1  =  Y,  7“so  ®  mo ,  (3-2) 

a. 

with  7a,  Sq  and  itlq  the  shearing  rate,  initial  slip  direction  and  initial  slip  plane 
normal  for  glide  system  a.  The  latter  two  vectors  deform  with  the  lattice  according 
to 

s“  =  FEsgS  m“  =  FE~Tm%,  (3.3) 

and  since  Sq  •  m]]  =  s“  •  m“  =  0,  plastic  deformation  is  isochoric 

jp  =  Jptr(FpFp_1)  =0,  Jp  =  1,  J=  JE.  (3.4) 

Ignored  are  possible  contributions  to  residual  deformation  and  plastic  volume  change 
from  nonlinear  elastic  and  core  effects  that  may  be  non-negligible  when  defect  den¬ 
sities  become  very  large  [Teodosiu,  1982;  Clayton  et  al,  2014].  Such  effects  could  be 
incorporated  explicitly  via  insertion  of  a  third  term  in  decomposition  (3.1)  [Clay¬ 
ton  et  al. ,  2005;  Clayton  et  al.,  2014],  but  are  not  considered  further  in  this  work, 
which  relies  on  the  traditional  two-term  kinematic  decomposition  of  crystal  plastic¬ 
ity  [Asaro,  1983].  Lagrangian  and  Eulerian  strain  tensors  of  (2.3)  are  redefined  here 
as  their  purely  elastic  counterparts 

E  =  i(FBTFp  -  1),  F>  =  ^(1-Fe-1Fe~t).  (3.5) 

The  usual  continuum  balance  laws  in  (2.4)  still  apply,  again  for  the  adiabatic 
case.  Respective  thermodynamic  potentials  for  Lagrangian  and  Eulerian  theories 
are  now 

U  =  U(E,  r],  C),  U  =  U{ D.ij.C);  tt  =  *(E,0,C),  tt  =  tt(D,0,C).  (3.6) 

Scalar  internal  state  variable  £  accounts  for  energy  of  lattice  defects,  e.g.,  stored 
energy  of  cold  work  associated  with  microscopic  stress  fields  of  dislocations  and  their 
core  energies  [Rosakis  et  al.,  2000;  Clayton  et  al.,  2014].  Replacement  of  (  with  a 
vector  or  higher-order  tensor  state  variable(s)  poses  no  conceptual  difficulties,  as 
has  been  demonstrated  elsewhere  for  Lagrangian  crystal  elastoplasticity  [Luscher 
et  al.,  2013]. 
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For  Lagrangian  elasticity  theory,  thermoelastic  constitutive  relations  for  stress 
(3.7)i,  temperature  (3.7)2;  entropy  (3.7)3;  conjugate  force  to  defect  density  (3.7)4; 
the  temperature  rate  form  of  the  local  energy  balance  (3. 8)1;  the  reduced  entropy 
inequality  (3.8)2;  and  the  conjugate  driving  force  (resolved  Kirchhoff  stress  [Asaro, 
1983])  for  slip  (3.9)  are  derived  from  (2.4)  and  (3.1)-(3.6)  using  standard  principles 
[Rosakis  et  al.7  2000;  Clayton,  2011;  Luscher  et  al. ,  2013] 

dU  d 'F  —*-1  -p-t  d  =  du/drij  r,=  -d^/d9:  x=-dV/d(, 


S  = 


dE 


I  - 


as 


c0  =  J2rar  +  0^:E 


(3.7) 

*-*!)* 

5>“t“  +  xC>o, 

a 

(3.8) 

=  (FetFeS) 

:  Sq  ®  mg  . 

(3.9) 

de 

a 

f"  =  Jer  :  sc 

For  Eulerian  elasticity  theory,  the  analogous  thermodynamic  relationships  are 
derived  as 

S  =  ^  =  jfet*fe7  e  =  dU/dri,  V  = -8^/86,  x  =  -a'F/ac, 


f)  c 

c0  =  ^f“7“  +  0-:D 


X 


(3.10) 

-0§f)c;  Ef^a  +  ^>°,  (3-11) 

'  rv 


ra  =  Jcr  :  sa  (g)  ma  =  (SF^  1FE-T)  :  (g  m£.  (3.12) 

On  the  right-hand  sides  of  energy  balances  (3. 8)1,  and  (3.11)i,  the  first  terms  rep¬ 
resent  plastic  dissipation,  the  second  terms  represent  thermoelastic  coupling  due 
to  thermal  expansion/contraction,  and  the  rightmost  terms  represent  stored  energy 
accumulation  with  defects  generated  during  deformation,  for  example.  As  defined 
in  (3.7)4  and  (3.10)4,  x  is  the  conjugate  thermodynamic  force  to  internal  state  vari¬ 
able  C  [Rice,  1971].  The  product  of  x  and  the  time  rate  of  £  gives  the  free  energy 
dissipation  rate  (if  positive)  or  free  energy  storage  rate  (if  negative)  associated  with 
evolution  of  the  local  defect  arrangement  (such  as  dislocation  density)  and  corre¬ 
sponding  residual  stress  state. 

Governing  equations  for  the  two  elasticity  representations  are  remarkably  simi¬ 
lar,  the  only  crucial  differences  being  the  different  stress  tensors  in  hyperelastic  laws 
of  (3.7)i  and  (3.10)i  and  resulting  differences  in  thermodynamic  driving  forces  for 
slip  in  (3.9)  and  (3.12).  It  is  remarked  that  quantitative  differences  in  resolved  shear 
stresses  (3.9)  and  (3.12)  arise  solely  from  truncation  of  Taylor  polynomials  corre¬ 
sponding  to  (3.6)  at  a  given  order  in  strain  measure  E  or  D.  For  elastic  strains  that 
are  not  too  large  (e.g.,  magnitude  less  than  unity),  such  Taylor  polynomials  should 
converge  to  give  the  same  values  of  resolved  shear  stress  if  a  sufficient  number  of 
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terms  (i.e.,  a  sufficient  number  of  higher-order  elastic  constants)  is  included  in  each 
polynomial.  Physically,  the  definition  of  each  of  (3.9)  and  (3.12)  is  identical:  The 
shear  stress  (current  force  per  unit  reference  area  or  energy  per  unit  reference  vol¬ 
ume  [Asaro,  1983])  acting  in  the  direction  of  slip  (or  parallel  to  the  Burgers  vector) 
and  resolved  on  the  plane  of  dislocation  glide. 

General  theories  are  complete  upon  prescription  of  kinetic  equations  for  inelastic 
deformation  and  defects,  e.g.,  for  Lagrangian  theory 


ja  =  7“(f“;  9,  C), 

C  =  C(x;E  ,9). 

(3.13) 

For  Eulerian  elasticity, 

r=r(fa;9,o, 

C  =  C(x;D  ,9). 

(3.14) 

The  first  of  (3.13)  and  (3.14)  apply  for  viscoplastic  slip.  For  idealized  rate  inde¬ 
pendent  plasticity  with  slip  resistance  (shear  strength)  ga  depending  possibly  on 
temperature  and  defect  content,  these  slip  rate  equations  could  be  replaced  with 
[Asaro,  1983] 

f“  <«7“(0,C)^7“  =  O,  f“  =  3“(0,C)^|7“|>O,  (3.15) 

r“  <ga(0X)^r  =  O,  Ta  =  ga(0,O^\ja\>0-  (3-16) 


Expanding  internal  energy  functions  as  Taylor  polynomials  of  order  three  in 
respective  elastic  strain  measures,  order  two  in  entropy  and  £,  and  bilinear  in  the 
product  of  strain  and  entropy,  (3.6)  produces 


U  —  ^  T  ^-'a(3x-^otEfiEx  T  OqT] 

U  —  —QapDaD p  -(-  ~C a/3^DaD T  Oqt/ 


2c 


+  2K  ’ 


1  -  Tn  Dn 


2c 


+ 


where  summation  applies  over  Greek  indices  denoting  Voigt  notation 
1,2,...,  6).  Second-  and  third-order  isentropic  elastic  constants  are 


(3.17) 

(3.18) 
{a  = 


d2U 

d2U 

r  d3U 

8EadE0 

o  dDadD0 

’  a0x  dEadE3dEx 

(3.19) 


^  d3U 

a0x  ~  dDadDgdDx  ’ 

with  a  zero  subscript  denoting  evaluation  at  the  unstressed  reference  state  at  ambi¬ 
ent  temperature  9q.  Relationships  between  third-order  Lagrangian  and  Eulerian 
constants  in  (3.19)2  and  (3.19)3  are  analogous  to  those  of  the  purely  elastic  the¬ 
ory  derived  in  (A. 8)  of  the  Appendix.  Ambient  specific  heat  at  fixed  strain  and 
Griineisen  parameters  are  equivalent  in  Lagrangian  and  Eulerian  representations 

c  =  dU/d6\0  =  dU/dO\0,  90Ta=  -d2U/dEadr]\o=  -d2(j/dDadg  |„.  (3.20) 

Residual  energy  of  lattice  defects  is  represented  by  the  rightmost  terms  in  (3.17) 
and  (3.18).  When  (  is  dimensionless,  /i  can  be  taken  proportional  to  a  representative 
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shear  modulus  Go-  For  example,  letting  if  denote  the  total  dislocation  line  length 
per  unit  volume,  b  the  magnitude  of  a  typical  Burgers  vector,  and  k  a  dimensionless 
parameter  of  typical  order  unity,  a  physically  realistic  assumption  is  [Clayton,  2011] 

C  =  byftj),  =  ^«G,  Qb2(f).  (3.21) 

3.2.  Nonlinear  elastic  effects  on  yielding  and  pressure  in  shock 
compression 

For  crystals  of  symmetric  orientations'3  subjected  to  planar  shock  loading  above 
the  HEL,  a  two- wave  structure  appears  consisting  of  an  elastic  precursor  followed 
by  a  slower-moving  plastic  wave  (presuming  the  shock  is  not  overdriven),  both 
waves  characterized  by  purely  longitudinal  particle  motion.  The  shock  stress  at 
the  HEL,  according  to  the  model  of  Sec.  3.1  and  written  as  Ph,  is  attained  when 
slip  rates  in  (3.13)-(3.16)  first  become  nonzero.  In  rate  independent  theory,  this 
corresponds  to  resolved  shear  stresses  of  (3.9)  or  (3.12)  attaining  critical  values 
ga.  Assuming  such  critical  strength  values  are  known,  the  model  of  Sec.  3.1  can 
be  used  to  analytically  determine  Ph,  and  the  corresponding  volume  ratio  at  the 
HEL,  Vh/Vo,  or  vice-versa,  assuming  no  plastic  deformation  or  defect  generation 
takes  place  for  (thermo) elastic  loading  up  to  the  HEL.  Similar  nonlinear  elastic 
analyses  have  been  used  to  infer  material  strength  from  experimental  longitudinal 
shock  data  [Lang  and  Gupta,  2010;  Turneare  and  Gupta,  2011].  The  importance  of 
using  nonlinear  elastic  rather  than  linear  elastic  models  for  such  calculations  has 
been  demonstrated  elsewhere  [Turneare  and  Gupta,  2011].  In  what  follows  next, 
differences  in  strength  predictions  resulting  from  different  choices  of  nonlinear  elastic 
model  (i.e. ,  Lagrangian  versus  Eulerian)  are  computed  for  several  representative 
cases. 

In  order  to  restrict  the  number  of  parameters  entering  the  analysis,  an  ideal 
solid  with  isotropic  symmetry  in  second-order  elastic  constants  and  simultaneous 
isotropic  and  Cauchy  symmetry  in  third-order  elastic  constants  is  now  considered 
[Clayton,  2013].  Cauchy  symmetry  of  third-order  constants  is  representative  of  sev¬ 
eral  noble  metals  [Hiki  et  al. ,  1966].  Isentropic  second-order  elastic  constants  are, 
in  terms  of  bulk  modulus  Bo  and  Poisson’s  ratio  v 

Cn  =  3— — Bo,  C12  =  3— — B0,  C44  =  Go  = - B0.  (3.22) 

1  +  v  1  +  v  2  1  +  v 

Isentropic  third-order  constants  are,  in  terms  of  Bo  and  its  pressure  derivative  Bq 

27  „  „  27 

Cm  =  5Cn2  =  I5C123  =  — ^-BqBo,  Cm  =  5Cn2  =  I5C123  =  — ^-(Bq  —  4)  Bo- 

(3.23) 

bFor  orientations  of  lower  symmetry,  the  response  to  step  loading  in  normal  stress  generally  consists 
of  quasi-longitudinal  elastic  and  plastic  waves  and  multiple  quasi-transverse  elastic  and  plastic 
waves  [Johnson,  1972]. 
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Uniaxial  strain  conditions  for  plane  wave  propagation  in  the  Ai-direction  are 
again  considered,  where  up  to  the  HEL,  with  tpo  denoting  the  initial  dislocation 
density 

Fp  =  1,  J  =  =  Fh  =  Vy  Vo,  (p  =  ipo  =  constant.  (3.24) 

Equations  (2.15)-(2.23)  still  apply,  but  many  reduce  to  simpler  forms.  Internal 
energy  functions  (3.17)  and  (3.18)  become 


U=  - 


3  1  -v 


--B'0E  )B0E2-d0(TE-l)V+- 


31  —  2v 


U  = 


21  +  v  14  0 


-- - -  -  —  (Bq  -4)£> 

21 +  v  14 v  0  ; 


k52B0+o, 


B0D2-  0o(TD-l)r]  + 


4  l  +  i/ 

31  —  2v 


(3.25) 


4  l  +  i/ 


Kb2B0ip0,  (3.26) 


where  the  quadratic  term  in  entropy  is  omitted  for  purely  elastic  shocks  [Thurston, 
1974], 

To  order  three  in  strains,  longitudinal  thermodynamic  stress  components  are 


Sir  =  3 


,1-"  27  D, 


Sii  = 


1  +  v  14 
1-1/  27 


-  —B'0E  1  Bn  E 


9 

28  ‘ 


Bn 


3  1  —  u 
41  +  1/ 


TBnE3, 


1  +  v  14 


-T7(Bo-4)^ 


Bn  D 


9  91-i/ 

—  (Bn  -  4)  H - 

28 V  0  ;  41  +  i/ 


(3.27) 


TB0D3.  (3.28) 


Terms  cubic  in  either  strain  measure  result  from  entropy  production  across  the 
shock.  Nonzero  Cauchy  stress  components  are  then,  for  Lagrangian  theory 

l-i/  27 

»  D' 


crii  =  -P  =  {1+2E) 


1/2 


—  —BnE  Bn  E 


1  +  v  14 


,  f  JLr'  _  31~p 

.  28  B°  41  +  v 


TBnE3 


(3.29) 


-1/2 


&22  —  C  33  —  (1  +  2E) 


V  31—^ 

28  0  41  +  1/ 


l  +  i/  70 


TBnE3 


-  i^B'E  )  BnE 


(3.30) 


Analogously,  for  Eulerian  theory, 

ern  =  -P=  (1-2 E)3/2  { 


Zl-^^{B’n-4)D 


[[1  +  v  14' 


Bn  D 


9  .  9  1-i/ 

28(B°“4)  +  4lT^ 


TBnE3 


(3.31) 


& 22  =  CT33  =  (1  -  2E)1/2 

9 


3^-S(go-4  )£> 


(Bn  —  4)  + 

28 '  0  ;  41  +  v 


1  +  v  70 v  0 
9  1  ~~  ^  TBnD3 


Bn  D 


(3.32) 
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Resolved  shear  stresses  from  (3.9)  and  (3.12)  are,  respectively 
fa  =  (1  +  2  E) 


3TT7  -  i  B“E 


+ 


+ 


9_ 

28  ‘ 


Bn 


3  1  -V 

41+7 


TBnE3 


3t7  -  f  B»E  1  B»E 


i  ,  9  ,  3  1  —  y 

(  28B°  41  +  z/ 


tb^e3  {s^m^  +  S3TO3), 


(3.33) 


Ta  =  (1  -  2D) 
9 


+ 


3TT7-n^-4^ 

9  1-1/ 


BnZ? 


(B'  -  4)  + 

28 '  0  ;  41  +  i/ 


TB0ir  }  s?m' 


51  "61 


3 — — - — (Bq  —  4)D 

l  +  i/  7(T  0  ; 


Bn  D 


+ 


28 


(Bq  —  4)  + 


9  1  -V 


41  +  1/ 


rBnD3  !>  (s^TO^  +  S3 7713 ), 


noting  that  under  the  present  uniaxial  strain  loading  conditions 


Jcr  :  s° 


m“  =  Jcr  :  Sq  (g)  mj  =  J(<7ns“m“  +  CT22S2  m%  +  ). 


(3.34) 

(3.35) 


Temperatures  in  respective  Lagrangian  and  Eulerian  solutions  are  simply 

9  =  dU/dr]  =  d0{l-TE),  9  =  dU/dr]  =  0O(1  —  TD).  (3.36) 

Two  representative  slip  configurations  are  considered:  (i)  Octahedral  slip  in  a 
cubic  crystal  shocked  along  the  cube  axis  [100]  and  (ii)  prismatic  slip  in  a  hexagonal 
crystal  shocked  along  the  a-axis  [1210].  Both  of  these  configurations  represent  pure 
mode  directions,  whereby  symmetry  would  lead  to  purely  longitudinal  elastic  and 
plastic  waves.  For  initiation  of  octahedral  slip,  eight  of  the  twelve  {111}  (110)  type 
systems  in  a  face-centered-cubic  (FCC)  crystal  would  be  subjected  to  an  equal 
resolved  shear  stress  r“,  the  remaining  four  unstressed  [Johnson  et  ai,  1970].  The 
same  solution  would  apply  for  { 1 10} (111)  slip  in  a  body-centered-cubic  (BCC) 
crystal,  with  s"  and  m“  simply  interchanged  in  (3.33)-(3.35).  Representative  slip 
vectors  are  taken  as 


s0  =  (1/72,  -1/72,0),  m0  =  (1/73,1/73,1/73).  (3.37) 

For  initiation  of  prismatic  slip,  two  of  the  three  {1010}  (1210)  systems  would  be 
subjected  to  an  equal  resolved  shear  stress,  the  third  unstressed  [Johnson,  1974], 
and 


s0  =  (—1/2,73/2,0),  m0  =  (73/2,1/2,0).  (3.38) 
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Resolved  shear  stress  magnitudes  obtained  from  Lagrangian  theory  by  (3.33) 
and  Eulerian  theory  by  (3.34)  are  shown  in  Fig.  5  for  octahedral  slip  system  (3.37) 
and  in  Fig.  6  for  prismatic  system  (3.38).  In  each  case,  the  pressure  derivative  of 
the  bulk  modulus  and  the  Griineisen  parameter  are  fixed  at  representative  values 


VIV o 


Fig.  5.  Resolved  shear  stress  on  {111}(110)  systems  for  shock  loading  along  [100]. 


Fig.  6.  Resolved  shear  stress  on  (1010}(1210)  systems  for  shock  loading  along  [1210]. 
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of  4  and  2,  respectively,  and  Poisson’s  ratio  is  varied  from  0.1  to  0.4.  Shear  stresses 
are  normalized  by  the  elastic  shear  modulus.  Differences  in  predictions  of  Eule- 
rian  and  Lagrangian  theories  emerge  at  volumetric  compression  ratio  {V /Vo)  less 
than  0.99,  i.e.,  at  greater  than  about  1%  uniaxial  compression.  For  a  fixed  value 
of  Poisson’s  ratio  jz,  shear  stresses  predicted  by  Eulerian  elasticity  tend  to  exceed 
those  predicted  by  Lagrangian  elasticity.  Normalized  shear  stresses  tend  to  increase 
with  increasing  v  because  shear  modulus  Go  decreases  rapidly  relative  to  the  bulk 
modulus  as  Poisson’s  ratio  increases. 

Table  3  shows  resolved  shear  stresses  predicted  by  each  theory  at  the  same  value 
of  HEL  stress  Pr  ■  These  are  the  slip  strengths  that  would  be  extracted  from  plate 
impact  data  on  Pr.  For  HEL  stresses  on  the  order  of  1%  of  the  ambient  bulk  mod¬ 
ulus  or  less,  differences  between  Lagrangian  and  Eulerian  predictions  are  negligible, 
as  inferred  from  the  first  row  of  Table  3.  This  would  be  the  case  for  very  ductile 
metals,  for  example.  However,  differences  between  Lagrangian  and  Eulerian  theories 
become  significant  for  strong  solids  with  larger  elastic  limits.  The  case  Pr/Bo  =  0.1 
is  physically  representative  of  prismatic  slip  in  shock-compressed  sapphire  [Gra¬ 
ham  and  Brooks,  1971;  Clayton,  2009]  or  silicon  carbide  [Clayton,  2010].  The  case 
Ph/Bo  =  0.25  is  physically  representative  of  octahedral  slip  in  diamond  [Lang  and 
Gupta,  2010].  In  the  latter  case,  the  difference  in  predictions  of  the  two  elasticity 
models  is  on  the  order  of  20%.  Slip  resistance  also  often  depends  on  temperature. 
Temperature  rise  predicted  by  Lagrangian  and  Eulerian  theories  through  (3.36)  is 
compared  in  Fig.  7.  Discernible  differences  emerge  at  compressions  exceeding  3%. 
Cleavage  fracture  is  known  to  occur  on  preferred  planes  in  crystals,  and  perhaps 
in  preferential  directions  due  to  lattice  trapping  [Riedle  et  al,  1996].  Results  of 
foregoing  analysis  may  apply  to  such  fractures. 

For  solids  with  relatively  low  shear  strength  (e.g.,  ga  <  O.OlBo),  shock  stress  P 
above  the  HEL  does  not  greatly  exceed  Cauchy  pressure  p,  and  deviatoric  stresses 
become  negligible  with  respect  to  pressure.  This  is  the  case  for  many  pure  ductile 
polycrystalline  metals,  and  is  approximately  true  for  ductile  single  crystals  when 
stresses  are  sufficient  to  enable  slip  on  primary  and  secondary  systems  so  that  the 
response  is  reasonably  isotropic  [Johnson,  1974].  The  same  statements  would  apply 
to  brittle  solids  with  low  fracture  strengths  on  a  sufficient  number  of  planes.  A 
standard  assumption  for  isotropic  solids  under  such  conditions  is  that  the  shock 


Table  3.  Resolved  shear  stresses  at  HEL  stress  Pr  with  v  =  0.25,  B],  =  4  and  T  =  2. 


Ph/Bo 

Ta/ Go  octahedral 

ra/ Go  octahedral 

Ta/ Go  prismatic 

ra/Go  prismatic 

0.010 

0.0045 

0.0045 

0.0048 

0.0048 

0.025 

0.0112 

0.0114 

0.0118 

0.0120 

0.050 

0.0219 

0.0227 

0.0233 

0.0241 

0.100 

0.0423 

0.0451 

0.0449 

0.0478 

0.250 

0.0939 

0.1107 

0.0996 

0.1174 
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Fig.  7.  Temperature  rise  versus  volume  ratio  for  elastic  shock  loading. 


pressure  can  be  adequately  represented  by  a  scalar  EOS,  i.e.,  P  ~  p  [Germain 
and  Lee,  1973;  Jeanloz,  1989].  If  isentropic  conditions  are  assumed,  as  in  (2.8)  and 
(2.9),  then  differences  between  shock  pressure  P  and  isentropic  hydrostat  p  are  due 
to  shear  strength  and  entropy  production  (e.g.,  from  thermoelasticity  and  viscous 


1.00  0.95  0.90  0.85  0.80  0.75  0.70 

VIV o 

Fig.  8.  Pressure  under  hydrostatic  compression,  generic  isotropic  solid. 
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dissipation)  that  would  occur  during  the  shock  process.  Predictions  of  third-order 
Lagrangian  EOS  (2.8)  and  third-order  Eulerian  EOS  (2.9)  are  compared  in  Fig.  8 
over  a  physically  realistic  range  of  Bq  .  Differences  between  the  two  become  apparent 
at  compressions  exceeding  about  5%,  with  Eulerian  theory  producing  relatively 
higher  pressure  than  Lagrangian  theory  for  a  material  compressed  to  the  same 
specific  volume. 

Polycrystalline  aluminum  and  copper  have  dynamic  deviatoric  strengths  ga  < 
O.OlBo.  Comparison  of  shock  data  with  hydrostatic  compression  data  [Mao  et  al. , 
1978;  Greene  et  al. ,  1994]  demonstrates  that  a  scalar  EOS  should  suffice  for  describ¬ 
ing  the  response  of  these  metals  shocked  compressed  over  0.95  >  V/Vq  >  0.75.  Mag¬ 
nesium  exhibits  relatively  low  resistance  to  slip  and/or  twinning  in  various  directions 
[Clayton  and  Knap,  2011,  2013;  Wu  et  al. ,  2012]  and  thus  a  scalar  EOS  should  also 
be  appropriate  for  untextured  Mg  polycrystals  at  moderate  compressions.  Predic¬ 
tions  of  each  third-order  EOS  in  (2.8)  and  (2.9)  are  compared  with  experimental 
shock  compression  data  [Marsh  et  al,  1980]  for  Al  in  Fig.  9,  Cu  in  Fig.  10,  and 
Mg  in  Fig.  11.  Compressibility  properties  of  Table  4  used  in  these  predictions  are 
obtained  from  ultrasonic  experiments  at  small  compression  [Guinan  and  Steinberg, 
1974]  and  are  not  fit  to  the  shock  compression  data.  In  each  case,  superior  agree¬ 
ment  of  Eulerian  theory  with  the  data  is  apparent.  Jeanloz  [1989]  suggested  that 
an  Eulerian  EOS  should  be  more  accurate  than  a  Lagrangian  EOS  of  the  same 
order  for  most  pure  substances,  in  agreement  with  the  particular  metals  analyzed 
here. 


Fig.  9.  Predictions  of  Eulerian  and  Lagrangian  EOS  and  shock  compression  data  for  aluminum 
[Marsh,  1980]. 
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Fig.  10.  Predictions  of  Eulerian  and  Lagrangian  EOS  and  shock  compression  data  for  copper 
[Marsh,  1980]. 


Fig.  11.  Predictions  of  Eulerian  and  Lagrangian  EOS  and  shock  data  for  magnesium  [Marsh, 
1980]. 
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Table  4.  Isentropic  bulk  modulus 
(GPa)  and  its  pressure  derivative 
[Guinan  and  Steinberg,  1974], 


Material 

Bo 

Bo 

A1 

76 

4.42 

Cu 

137 

5.48 

Mg 

35.4 

3.90 

4.  Discussion 

Results  and  analysis  in  Secs.  2.3  and  3.2  suggest  that  overall,  Eulerian  theory 
should  provide  a  more  accurate  depiction  than  Lagrangian  theory  of  the  mechan¬ 
ical  response  of  metallic  crystals  under  spherical  compression  and  uniaxial  strain 
compression,  in  particular  when  internal  energy  functions  are  truncated  at  orders 
two  or  three  in  corresponding  finite  strain  measures.  This  assertion  can  be  under¬ 
stood  to  follow  from  the  physical  observation  that  in  typical  metallic  crystals,  strain 
energy  density,  stress  or  pressure,  longitudinal  stiffness,  and  tangent  bulk  modulus 
all  tend  to  increase  rapidly  with  decreasing  specific  volume,  and  from  the  math¬ 
ematical  observation  that  components  of  Eulerian  strain  increase  more  rapidly  in 
magnitude  in  compression  (become  more  largely  negative)  than  do  components  of 
Lagrangian  strain,  as  demonstrated  explicitly  in  Table  5.  For  example,  second-order 
Eulerian  elasticity  tends  to  predict  significantly  larger  (in  magnitude)  axial  stress 
than  second-order  Lagrangian  theory  since  in  second-order  theories,  thermodynamic 
stresses  are  linearly  related  to  work  conjugate  strains.  It  is  clarified  that  the  results 
and  analysis  herein  do  not  necessarily  favor  one  strain  measure  over  the  other  for 
general  descriptions  of  continuum  kinematics.  Rather,  they  support  the  notion  that 
a  finite  Taylor  series  representation  of  internal  energy  in  Eulerian  tensor  D  appears 
advantageous  to  one  of  the  same  order  in  Lagrangian  tensor  E  for  describing  the 
stress  response  of  possibly  anisotropic  metallic  materials  under  finite  compressive 
strain.  Furthermore,  although  such  Taylor  series  representations  of  energy  are  math¬ 
ematically  rigorous  and  convenient  for  general  anisotropic  solids  since  they  incorpo¬ 
rate  correct  symmetries  through  corresponding  independent  elastic  constants,  such 
polynomial  forms  may  not  necessarily  be  ideal  for  describing  particular  materials. 
Another  drawback  for  such  higher-order  theories  is  that  elastic  constants  of  order 


Table  5.  Axial  strain  components  under 
spherical  and  uniaxial  deformation. 


V/V0 

Spherical  strain 

Uniaxial  strain 

Eu 

Du 

Eu 

Du 

0.7 

-0.106 

-0.134 

-0.255 

-0.520 

0.8 

-0.069 

-0.080 

-0.180 

-0.281 

0.9 

-0.034 

-0.036 

-0.095 

-0.117 

1.0 

0.000 

0.000 

0.000 

0.000 

1.1 

+0.033 

+0.031 

+0.105 

+0.087 
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three  are  difficult  to  measure  very  accurately  in  experiments  [Thurston,  1974],  and 
those  of  orders  four  and  higher  are  even  more  difficult  to  obtain  and  have  been 
reported  for  very  few  crystalline  solids. 

While  such  benefits  of  Eulerian  theory  are  readily  apparent  for  loading  condi¬ 
tions  relevant  to  shock  compression  and  hydrostatic  compression,  the  same  asser¬ 
tions  do  not  necessarily  hold  for  other  loading  protocols.  In  tension,  for  example, 
Lagrangian  strain  components  increase  more  rapidly  than  Eulerian  components 
(Table  5),  opposite  to  what  occurs  in  compression,  leading  to  different  trends  in 
predicted  material  response.  However,  most  crystalline  metallic  samples  of  signifi¬ 
cant  size  will  undergo  cavitation  or  tensile  fracture  when  subjected  to  finite  tensile 
volumetric  strain,  and  they  will  undergo  plastic  slip,  twinning,  or  shear  fracture 
when  subjected  to  large  deviatoric  strain.  Thus,  effects  of  choice  of  Lagrangian  ver¬ 
sus  Eulerian  theory  become  of  low  significance  since  large  elastic  tensile  or  deviatoric 
deformations  are  not  sustainable.  On  the  other  hand,  in  other  material  classes  such 
as  rubber  and  ductile  polymers,  large  elastic  tensile  and/or  shear  deformations  are 
often  sustainable,  and  nonlinear  elastic  potentials  (neo-Hookean,  Mooney- Rivlin, 
etc.)  that  may  differ  significantly  from  the  present  Lagrangian  or  Eulerian  Taylor- 
type  polynomials  have  been  successful  for  describing  their  responses,  often  isotropic 
and/or  (nearly)  incompressible  [Ogden,  1984].  Large  elastic  deviatoric  strains  may 
also  be  possible  in  some  molecular  crystals  [Clayton  and  Becker,  2012]  or  in  other 
crystalline  nonmetals  with  restricted  slip  systems,  especially  under  conditions  where 
confining  pressures  are  large  enough  to  suppress  fracture,  for  example  indentation. 

5.  Conclusions 

The  response  of  metallic  solids  to  planar  shock  loading  has  been  analyzed  in  the 
context  of  nonlinear  Lagrangian  and  Eulerian  thermoelasticity  models,  the  latter 
incorporating  a  relatively  uncommon  “Eulerian”  strain  referred  to  material  coordi¬ 
nates.  Predictions  of  shock  stress  for  Eulerian  and  Lagrangian  theories  have  been 
compared  for  several  anisotropic  metallic  single  crystals  (Al,  Cu  and  Mg),  with  elas¬ 
tic  constants  of  up  to  order  four  included.  Differences  between  predictions  of  each 
model  as  the  number  of  higher-order  constants  is  increased  are  smaller  for  Eulerian 
theory  than  Lagrangian  theory,  and  higher-order  constants  themselves  are  of  smaller 
magnitude  for  Eulerian  than  Lagrangian  theory,  suggesting  greater  accuracy  and 
faster  convergence  of  Eulerian  over  Lagrangian  theory.  Eulerian  theory  also  more 
closely  matches  atomic  simulation  data  for  elastic  shock  compression  of  aluminum. 

A  continuum  crystal  plasticity  framework  has  been  developed  that  incorporates 
an  Eulerian  elastic  strain  tensor  in  the  thermodynamic  potentials.  The  fundamental 
difference  between  this  framework  and  traditional  Lagrangian  crystal  elastoplastic- 
ity  in  the  resulting  thermodynamic  relations  is  the  different  constitutive  equality  for 
stress,  which  can  result  in  different  driving  forces  for  slip  (i.e. ,  critical  resolved  shear 
stresses)  and  different  pressure-volume  behavior  when  internal  energy  potentials  are 
truncated  at  finite  order.  Differences  in  critical  resolved  shear  stresses  derived  for 
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the  two  third-order  elastic  models  have  been  shown  to  become  important  at  vol¬ 
umetric  compressions  in  excess  of  about  one  percent.  Different  values  of  resolved 
shear  stress  would  be  obtained  from  the  two  different  models  from  axial  shock  Hugo- 
niot  data  for  materials  with  Hugoniot  elastic  limits  in  excess  of  several  percent  of 
their  bulk  modulus,  an  important  consideration  for  deducing  dynamic  strengths 
of  strong  solids.  At  shock  pressures  above  the  HEL  in  polycrystalline  Al,  Cu  and 
Mg,  Eulerian  theory  provides  a  better  fit  to  experimental  shock  compression  data 
than  Lagrangian  theory  of  the  same  order.  The  analysis  empirically  suggests  that 
this  Eulerian  theory  is  preferable  to  traditional  Lagrangian  theory  for  modeling  the 
shock  response  of  these  metals. 


Appendix:  Elastic  Potentials,  Elastic  Constants,  and 
Equations-of- State 

Consider  Taylor  polynomial  expansions  of  Lagrangian  and  Eulerian  internal  energy 
functions  in  (2.5)  of  order  three  in  respective  strain  measures  E  and  D,  restricted 
here  to  isentropic  conditions  (i.e.,  omitting  77  dependence)  for  illustrative  simplicity 

U(E)  =  -C  ijklEijEkl  +  -Cijklmn  Eu  EklEMn  +  0(||E||4),  (A.l) 

U(  D)  =  -C  ijklDijDkl  +  -CijklmnDjjDklDmn  +  0(|jD||4).  (A. 2) 

Terms  linear  in  either  strain  measure  vanish  by  presuming  that  the  unstrained 
reference  state  is  stress  free,  and  strain  energy  is  also  assumed  to  vanish  in  this 
unstrained  reference  state.  The  following  kinematic  identities  apply: 


dEu/dFkL  =  SilFw  +  SjLFki), 

dDu/dF \L  =  \f£{F£FJ^  +  FJ^F^). 

Chain  rule  differentiation  and  use  of  (2.6),  (2.7)  and  (A. 3)  then  lead  to 
d2U  „  „  d2U 


=  FiNFkM 


dEJNdELM 


—  p-l  p-l  p-l  p-1  p-1  p-1 

-  ±  1  TU  1  Tm  -L  D™  T  r  m 


d2U 


dFijdFkL 
d2U 

dFijdFkL  ~  ±  Qi  ±  Ik  Jm ±  Rm*  Ln  ±  Pn  QD qrQD Ip 

~  FJkP'L  ~  FUPkJ  -  P7™FLmFkNPzN, 


(A.3) 


+  ' hkFj^PnL,  {PiJ  =  JF jkCTik),  (A. 4) 


(A.5) 


with  P  the  first  Piola-Kirchhoff  stress.  In  the  undeformed,  stress  free  reference 
state,  (A. 4)  and  (A.5)  reduce  to 


/  d2U  \ 


(  d2ti  \ 


\dFijdFkL)  p=1  S'M5kNCMJNL ’  1  dF.udFkL  I 


F=1 


(A.6) 
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Assuming  that  second-order  tangent  moduli  (second  derivatives  of  internal  energy 
with  respect  to  F)  are  equal  at  the  reference  state  by  equating  left-hand  sides  of  each 
of  (A. 6),  it  follows  that  second-order  elastic  constants  are  equivalent  in  Lagrangian 
and  Eulerian  theories  (Voigt  indices  in  Greek) 

C IJKL  =  c  IJKL  =  c  IJKL,  Cc</3  =  Ca/3  =  Ca/3-  (A. 7) 

Extending  the  same  procedure  by  differentiating  once  further  with  respect  to  F  and 
equating  third-order  moduli  d3U/dFdFdF  for  Lagrangian  and  Eulerian  represen¬ 
tations  in  the  stress  free  reference  state,  third-order  elastic  constants  are  related  by 
[Clayton,  2013] 


C IJKLMN  =  ClJKLMN  +  SlK^JLMN  +  &IlS~JKMN 

+  Sim^kljn  +  SinCkljm  +  SjkCilmn 
+  SjmCinKL  +  SjlCikMN  +  SjnCimkl 
+  SkmCijln  +  SknC-ULM 

+  SlmCijkn  +  SlnC-ijkm-  (A.8) 

Extending  this  procedure  to  fourth  order  results  in  (2.25)  for  the  uniaxial  strain 
case.  The  second  of  (2.14)  is  derived  by  including  polynomial  terms  in  entropy 
change  in  internal  energy  functions  (A.l)  and  (A. 2)  and  equating  cross-derivatives 
of  internal  energies  with  respect  to  F  and  rj  [Clayton,  2013].  Note  that  (A. 7)  is 
consistent  with  the  requirement  that  when  strains  are  small  enough,  the  difference 
between  E  and  D  becomes  negligible  (both  approach  the  strain  tensor  of  linear 
elasticity),  and  the  second-order  elastic  constants  become  identical  to  those  used  in 
linear  elasticity.  Result  (A.8)  is  consistent  with  the  relationship  derived  elsewhere 
[Perrin  and  Delannoy,  1978]  by  substituting 

D=E-2E2+4E3 - «E-2E2  (A.9) 

into  (A. 2)  and  then  equating  terms  of  order  three  in  E  in  (A.l)  and  (A. 2).  Note  that 
if  the  Eulerian  strain  energy  potential  (A. 2)  is  exact  at  order  three,  for  example,  then 
the  difference  between  internal  energy  functions  (A.l)  and  (A. 2)  is  the  truncation 
error  of  the  former,  which  is  0(||E||4).  In  other  words,  imposition  of  (A. 7)  and 
(A.8)  does  not  ensure  that  (A.l)  and  (A. 2)  are  equal  since  terms  of  orders  four 
and  higher  in  strain  will  lead  to  differences  between  the  two  energy  functions  that 
tend  to  become  more  pronounced  at  larger  strain.  However,  if  (A.l)  and  (A. 2)  are 
extended  to  polynomials  of  order  n  in  strain,  then  truncation  errors  are  of  order 
n  +  1  in  some  norms  of  strain  tensors  E  and  D;  these  errors  tend  towards  zero  as 
n  — >  oo  for  regimes  in  which  such  norms  are  of  less  than  unit  magnitude,  in  which 
case  (A.l)  and  (A. 2)  become  both  exact  and  equivalent. 
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For  materials  of  high  cubic  or  isotropic  symmetry,  the  bulk  modulus  and  its  first 
derivative  with  respect  to  pressure,  in  the  reference  state,  are,  in  Voigt  notation 


Bo  —  g  (Cn  +  2C12), 


(A.10) 


(1  -  - 

2-  \ 

/  1  1 

(  1  ~  ~ 

2  -  \ 

/ 

f  -Cm  +  2C112  + 

3Cl23l 

/Bo  =  -  3  i 

,  3  Cm  +  2Cn2  + 

3Cl23l 

/Bo  +4, 

(A.ll) 


Derivations  of  the  first  of  (A.ll)  can  be  found  in  Thurston  [1965]  and  Teodosiu 
[1982].  The  second  equality  of  (A.ll)  then  follows  upon  use  of  (A. 8).  Letting  F  = 
J 1/31  in  (A.l)  and  (A. 2),  and  then  applying  stress-strain  relations  (2. 6)1  and  (2.7)i, 
Cauchy  stress  becomes  of  the  spherical  form  <x  =  —pi,  where  the  resulting  pressures 
p  corresponding  to  (A.l)  and  (A. 2)  are  given  by  equations-of-state  (2.8)  and  (2.9), 
respectively.  Truncation  errors,  which  are  not  written  explicitly  in  these  equations- 
of-state,  are  roughly  of  order  three  in  corresponding  strains  since  truncation  errors 
in  internal  energy  are  of  order  four. 
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